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Abstract. We introduce a two-state non-conserving driven-diffusive system in 
one-dimension under a discrete-time updating scheme. We show that the steady- 
state of the system can be obtained using a matrix product approach. On 
the other hand, the steady-state of the system can be expressed in terms of a 
linear superposition Bernoulli shock measures with random walk dynamics. The 
dynamics of a shock position is studied in detail. The spectrum of the transfer 
matrix and the relaxation times to the steady-state have also been studied in the 
large-system-size limit. 
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1. Introduction 

Exactly solvable systems have been of great interests to physicists for years. Among 
these systems there are one-dimensional non-equilibrium systems which have unique 
critical and collective properties [J- [I]- Some of these properties, which usually 
can not be found in their equilibrium counterparts, are out-of-equilibrium phase 
transitions and shock formations. Despite of these interesting properties, the number 
of one-dimensional driven-diffusive systems belonging to the family of non-equilibrium 
systems, which can be solved analytically is very limited. 

Over the past two or three decades, the one-dimensional driven-diffusive systems have 
been studied extensively from different point of views. It has been shown that some 
of these systems can be mapped onto the zero-range processes [3]. Some of them, on 
the other hand, are related to the lattice-path models [IJ [5l |6] . It is known that in 
some of these systems the dynamics of a shock (or equivalently a sharp discontinuity 
in the density of particles in the system) is similar to that of a simple random walker. 
In this way the study of a system with large number of degrees of freedom reduces to 
the study of a system with fewer number of degrees of freedom [7] . 
It has been shown that in the continuous-time updating scheme there are only 
three two-state driven-diffusive systems with open boundaries and nearest-neighbor 
interactions in which a single shock can evolve in time with a random walk dynamics 
[7]. The steady-state of these systems can be written as matrix product states and 
that their associated quadratic algebras have two-dimensional matrix representations 
[8]. This idea has been also generalized to the systems containing multiple species of 
particles with next-nearest-neighbor interactions [9l[ini[lli. However, in the discrete- 
time updating scheme there is no such a classification and only a couple of examples 
exist [lailS]. 

In this paper we introduce a two-states one-dimensional driven-diffusive system with 
open boundaries. The updating scheme is parallel which consists of two steps. This 
scheme is sometimes called the sublattice-parallel updating scheme. This system does 
not belong to that three-member family of two-states driven-diffusive systems with 
open boundaries and nearest-neighbor interactions. The steady-state of this system 
can be obtained using a matrix product approach (for a review see [14]). On the other 
hand, the steady-state of the system can be written as a superposition of shocks with 
simple random walk dynamics. The time evolution of the probability distribution 
function can also be studied in special cases. The eigenvalues and eigenvectors of 
the time evolution operator can be obtained using a simple plane wave ansatz. The 
relaxation times can be calculated exactly for large-system-size. 

This paper is organized as follows: in the first section we define the model. By 
introducing two product shock measures we show that the time evolution of the shock 
positions are similar to those of two simple one-dimensional random walkers in discrete 
time. We build the steady-state of the system as a linear superposition of these shocks. 
In the second section we review the basics of the matrix product approach and show 
that the steady-state of the system can also be found using this approach. In the 
third section we investigate the relaxation to the steady-state. In this direction we try 
to find the largest eigenvalue of the time evolution operator using the standard plane 
wave ansatz by imposing some restrictions on the microscopic reaction probabilities. 
We will finally discuss about the possible extensions to this model and open questions. 
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2. Definition of the process 

We consider a lattice of length 2L with open boundaries. Each lattice site can be either 
occupied by a particle ^ or a vacancy 0. The particles are injected and extracted from 
both ends of the lattice with certain probabilities. At left boundary (the first lattice 
site) we have: 

$ ~> A with probability a 

A with probability 7 
while at right boundary (the last lattice site): 

^ with probability S , , 

^ ^ with probability /3. ^ ' 

We assume that in the bulk of the lattice the particles hop to the left and right 
while interacting with each other. Apart from diffusion we have pair annihilation and 
coagulation of particles to the left and right. The reaction rules are as follows: 

+ A ^ A + with probability 1 

^ + ^ + A with probabihty 1 

yl + A ^ + with probability tu (3) 
A + A ^ + A with probabihty t24 
A + A ^ A + ib with probabihty t^^ 
under the constraint: 

tl4 + t2i + i34 = 1 (4) 

which will be discussed in the next section. The time is discrete and that the updating 
scheme applied is called the sublattice-parallel which is defined as follows: we divide 
the bulk dynamics into two half-time steps. In the first half-time step the pairs of 
neighboring sites {2k, 2k +1) for k = 1, ■ ■ ■ , L — 1 and also the first and the last 
lattice sites are updated. In the second half-time step the pairs of neighboring sites 
{2k — 1, 2k) for k = 1, • • • , L are updated. 

The time evolution of the probability distribution vector \P{t)) is given by the following 
master equation: 

T\P{t))^\P{t+l)). (5) 
where T in ([5]) is called the transfer matrix. The transfer matrix T is given by a 
product of two factors Ti and r2 as T = TiT2. These factors are defined as follows: 

Ti = c®T® ...®T®n = c®T®'^^~^'^ ®n 
T2^T(E)T(E)...(E)T = r®^ 

where T, C and TZ are given by: 
/I ii4 \ 



r 



). (6) 



a 1-77' V ^ 1-/3 



1 t24 
1 ^34 
\ / 

The matrix T is written in the basis (00, 0A, A0, AA). The matrices C and TZ are also 
written in the basis (0, A). In the long-time limit the system approaches its steady- 
state and one has: 

T\P*) = \P*). (7) 
In the following sections we will investigate the steady-state of the system using two 
different approaches. 
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3. Steady-state as a superposition of shocks 

In this section we sliow tliat one can construct ttie steady-state of tlie system as a linear 
superposition of product shock measures. At the outset we define two product shock 
measures and investigate their time evolution under the transfer matrix T defined in 
the previous section. The product shock measures at even sites 2k {k — 1, ■ ■ ■ , L) will 
be denoted by \2k) and the product shock measures at odd sites 2fc + 1 {k — 0, ■ ■ ■ , L) 
will be denoted by |2fc + 1). Defining the notation: 



o.e\ 



91.2 



\Pi:2) = 

\ Pl,2 

we write: 

\2k) = ® \pl}f''' ® ® IpDr"--" (8a) 

2k 

\2k + 1) = (|p?) ® ip^r'^ ®\p^^\p'2) ® ® iP2»^'^~'"'- m 

2k+l 

Note that the lattice site 2L + 1 is an auxiliary lattice site so that the product shock 
measure \2L + 1) indicates a distribution of particles with densities and pf at odd 
and even lattice sites respectively. In this case the shock front can be considered to 
be between the lattice sites 2L and 2L + 1. We now investigate the time evolutions 
of (|8ap and ((86| under the transfer matrix T. Let us consider the following values for 
the density of particles: 

pI^O , p°^6 , pl^a , pl = 0. (9) 

These values beside the constraint provide the necessary and sufficient conditions 
which guarantee that the shock position has a simple random walk dynamics. After 
some calculations we find: 

T|l) = ^,i|2) + ^,2|3) + (1 - TTrl - 7r,2)|l) 

T|2fc) =7r;|2fc- 1) +7r^|2fc + l) +7r,|2A;) for l<k<L 

T\2k + 1) = 7r;7r^|2fc) + 7rr7r,|2fc + 2) + TTl\2k + 3) + TTf\2k - 1) (10) 

+ {-Ks + 2iiiiir)\2k + 1) for l<k<L-\ 
T\2L + 1) = ^ii\2L) + 7Ti2\2L - 1) + (1 - nn - n2)\2L + 1) 

in which we have defined tTs = 1 — tt^ — tt; and that: 

TTr = 1 — (5 + (5^24 , TTi = 1 — a + a 

_ a(l-,5) + (l-7)^ _ a(l-^) + (l-7)g x 

TTrl — T^s 1 T^r2 — T^r U -L j 

S{l-a) + (l-l3)a S(l~a) + {l-B)a 

The equations (jlOl) indicate that the shock position evolves in time with a dynamics 
which is similar to that of a simple random walker and when it reaches to the 
boundaries of the lattice it refiects from them. 

Using the equations pUj) it is now easy to construct the steady-state of the system 
defined in JTJS]) . Let us consider a linear superposition of the above defined shocks as 
follows: 

, 2L+1 

m-^E^'cifc) (12) 
^ fc=i 
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in which Zj^ is a normalization factor given by: 

2L + 1 
k=l 



(13) 



If \P*) is the steady-state then it should satisfy ([7]). It turns out that the coefficients 
Pjf's are given by: 



p* 



^2 — 



2L 



^ TTi 

)) 

7rii7ri+-n"!2(l-7rs) ^ ttj • 



(14) 



PI* 
9 



' TTr \2L+1 



Note that if a 
can be written as: 



_ Hi fZtV 

2L+1 irii7ri+Ti2(l-7rs) "> tt, ^ 

(5 = then an empty lattice is the steady-state of the system which 



IP*) = 



®2L 



(15) 



On the other hand, ii a — (3 — 1 and < d,j < 1 the steady-state of the system is 
given by: 



(16) 



Finally, if 5 = 7 = 1 and < a, /? < 1 the steady-state of the system is given by: 



IP*) = 



l8i 



(17) 



In the next section we describe how one can obtain the steady-state of the system 
using a matrix product approach. 



4. Matrix product approach 

According to the matrix product approach the steady-state of the system |P*) can be 
written as [15) : 

-, ®L 



\p*) = ^^m\ 



E 

b 



E 
D 



\V)) 



(18) 



in which the operators E and D (E and D ) stand for the presence of a vacancy 
and a particle at an odd (even) lattice site respectively. The denominator Zl is a 
normalization factor. The four operators E, D, E, and D besides the vectors \V)) 
and ((W^l satisfy a quadratic algebra which results from the following relations: 



mc[ I ) = {{w\[ 



n 



E 
D 



)\V)) - ( 



E 

b 



E 
D 

\V)). 



(19) 
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The normalization factor is also given by: 

ZL = {{W\iCCf\V)) (20) 

in which C = D + E and C ^ D + E. Assuming C* ^ C by using dH) and (HH) we find 
the following quadratic algebra associated with the process (Il])-(|3]): 

DD = , [D, C] = (^34 - l)DD , [C, D] = {t^i - l)DD, 
{{W\{aC+il-a-j)D-D)^0, (21) 
((5C + (1 - /3 - S)D - D)\V)) = 0. 

The quadratic algebra (j2ip has a two-dimensional matrix representation provided that 
the constraints (U) and ^ are fulfilled. The representation is then given by : 





P°2 





d 


—p°l 


Pi 





d 





-d t(l-p?) 



1 -P2 



1 



^ +^ ^+ V 7r,(a(5+/3-l)-5) 



where d = — d( )■ Note that the algebra pip has also one-dimensional 
representations which are associated with the special cases discussed at the end of 
the previous section. In this case the operators E, D, E and D are in general complex 
numbers. For a = 5 = Q one finds D = D = and C = E ~ E. On the other hand, 
for a = /? = 1 and < 5, 7 < 1 one finds D = E = and C E = D. Finally, 
for (5 7 = 1 and < a, /3 < 1 one finds D ^ E ^ and C ^ D ^ E. Apart 
from these special cases, the two-dimensional matrix representation given above is the 
exact result. Since the steady-state of the system is unique, it is not difficult to show 
that the steady-state of the system in terms of a linear superposition of shocks given 
in (jl2p is exactly equal to one given in (jlSp in terms of a product of non-commuting 
operators. 



5. Relaxation to the steady-state 

In the following we study the relaxation time to stationarity in the random walk 
picture using a plane wave ansatz. Let us define Pk{t) as the probability of finding 
the shock position at time t at the lattice site k. As t oo the probability Pk{t) 
converges to given in p^ . We write: 

2L+1 

\P{t)) = E ^^Wl^)- (23) 

k=l 
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The time evolution of is governed by the master equation ([5]). Using (O, pU)) 
and and after some calculations we find: 

Pl(t + 1) - TTiP2it) + nfPait) + (1 - TTrl - 7r,2)Pi(i) 
P2{t + 1) = TTrlPlit) + TTlTTsPsit) + TT.Pa W 

Psit + 1) = ^,2^1 (i) + nrP2{t) + TTiPiit) + nfPsit) + (vr, + 2^i^,)P3(i) 
P2fe(t + 1) =7r,,7rsP2fc-i(i)+7ri7r^P2fe+i(t)+7r^P2fe(t) for 2<fc<L-l 

P2fe+l(t + 1) - TT^P2k-lit) + TTrP2k{t) + TTiP2k+2{t) + T^t P2k+z{t)+ ,^ 

(tt, +27r;7r,)P2fc+i(t) for 2<fc<i-2 ^^^> 

P2L-l{t + 1) = ^^P2L-3W + 7r,P2L-2W + nP2L{t) + TTl2P2L+l{t) + 

(vTs + 2'KlTTr)P2L-l{t) 
P2L{t + 1) = nrn,P2L-l{t) + nnP2L+l{t) + TlsP2L{t) 
P2L+l{t + 1) = 7r,P2L(0 + 7r2P2i_i(<) + (1 - TTa - 7ri2)P2L+l(i) 

Before going any farther let us have a note of the drift velocity v and diffusion 
coefficient D of the shock front in an infinite system. Considering the bulk equations 
in (I24p for an infinite system, one can calculate v and D according to the approach 
used in The moments are given by: 



oo 



by defining: 



E k^'Pk{t) = {K{t)) + {K{t)) (25) 

k— — oo 

(Kit)) ■■^EZ-o.m-P2kit) 

(Kit)) ■■^j:Z-ooi2k+irP2k+i{t). 

Using the definitions: 

v = \iint^oo{{k{t + l)) - {k{t)}) 

D = limt^U{{P{t + f )) - {kit + 1)>2) - ((fc2(t)) _ (fc(i))2)). 

one finds: 

^ = 2^ , D^2\-^il^f). (26) 

1 + TTs 1 + TTs 2 

As can be seen these expressions are exactly the same as those belong to the 
asymmetric simple exclusion process calculated in |12| . This means that in the long- 
time limit these properties do not depend on the microscopic dynamics. 
Finding a general solution for the equations (p4|) is a formidable task; therefor, in what 
follows we consider the case tt^ = which results in tTj-i — Tin — 0. In this case only 
the shocks at odd lattice sites appear in Simplifying the equations (1^^ results 

in: 

Pi(t + I) = Tr2p3(t) + (1 - 7r,2)Pi(0 

P3(t + I) = ^.2Pl(i) + 7rf P5(t) + 2TIi^rP3{t) 

P2k+iit + l) = nfP2k-i{t)+T:fP2k+3{t)+27rnrrP2k+i{t) for 2<fc<L-2 (27) 

P2L-l(i + 1) = 7TfP2L-3it) + ni2P2L+lit) + 2n^rP2L-l{t) 
P2L+l{t + 1) = 7r2P2L_i(<) + (I - 7r,2)P2L+l(t) 

where tt; + tt,. = 1. For the eigenfunction we suggest: 

Pk{z,t) ^ k\z)Pk{z) for odd k (28) 
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in which A(z) is the eigenvalue. By considering: 

P2k+i{z) - (f^)^(a(z)z2fc+i - a(z-i);j-2fe-i) for 1 < fc < L - 1 

P^{z) ^ {^)Ha{z)z - a{z-^)z-^) (29) 

P2L+i{z) = C2 (|^)^(a(z)z2^+i - 
and replacing them in (j27p and after some straightforward calculations one finds: 

a{z) = {^)kz~\T,,T,iz-^ + (1 - 7r,2 - A(z))Ci) , Ci = ^ , = ^ 

and the eigenvalues: 

A(z) =7r,.7r/(z + z-i)2. (30) 

It turns out that the proper z's belong to two groups: The z's of the first group 
contains four z^s which are z = ±(^)*2 resulting in A = 1 and therefore correspond 
to the steady-state. The z's of the second group satisfy the following equation: 
_^,,_ {1 + Az^){l + Bz^) 

(z2+^)(z2+B) ^ ' 

_ 2 _ 2 __ 

where A = '^'^^ J'' and B ~ '^'''^^J'' ■ The equation (PT|) has 4L + 4 roots; however, 
the roots z — ±1 and z = ±i result in a zero eigenfunction and have to be excluded. 
One also notes that since A(z) = A(— z) and also A(z) = A(i) the remaining 4L roots 
result in L eigenvalues. 

The roots of a very similar equation to ([31]) has already been investigated in |16| . 
Finding the roots of (I3ip for an arbitrary L can be quite difficult. However, for 
L — ^ oo the roots of pip are greatly simplified and the eigenvalues can be written as 
follows: 

A = 47r^7r, for - 1< AT, F < 3 + 2\/2 (32) 
A = -7r^7ri(Vy- -^f for Y > 3 + 2\/2 , -\ < X <Y (33) 




A = 7r^7ri(V|A:| + ^=)^ for X<-1, r>X (34) 

where X and F correspond to A and B ox B and ^. Now the largest relaxation time 

is given by t = |i?e[ln A]|~^. If A or B are such that in the thermodynamic limit one 

of the roots of (|3ip results in A = 1, then the relaxation time diverges. Note that the 

eigenvalues (|32t - (|34l) do not always occur and their presence depends on whether the 

hopping probabilities 7r;_r and the reflecting probabilities i^r2,i2 are positive. 

Let us consider a simple example for which the roots of (f3T|) can be calculated exactly. 

For: 

134 = J — , t24 = (o5) 

b a 

we have: 

n2 = 7rf , TTr2 = T^l ■ (36) 

In this case A ^ B ~ and the roots of ([31]) are, according to ((32|). phases and result 
in a relaxation time which in the large-system-size limit is given by : 

T = \\n{AnTTrT^. (37) 
At the phase transition point, where tt; = tt^ = ^ and the shock drift velocity is zero, 
the relaxation time is diffusive. 
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6. Concluding remarks 

In this paper we have introduced a family of driven diffusive systems on an open 
lattice under sublattice-parallel updating scheme. We have calculated the steady- 
state of the system using the matrix product approach with a two-dimensional matrix 
representation whose structure is quite similar to that introduced in [17] for the case 
when the steady-state of the system can be written in terms of a linear superposition 
of shocks with random walk dynamics. We have shown that the dynamics of a single 
shock in the system is similar to that of a single random walker which moves on a 
one-dimensional lattice with reflecting boundaries. We have calculated the hopping 
probabilities of the shock front in the bulk and also the probabilities of reflection from 
the boundaries. Using the random walk picture we have been able to calculate the 
relaxation times in the limit of large-system-size. 

The asymmetric simple exclusion process is the only model which has already been 
introduced and studied under sublattice-parallel updating scheme. Apart from the 
model studied in this paper defined in ([T])-([3]), we have also found two other processes 
containing only one species of particles under sublattice-parallel dynamics in which 
the random walker picture can be applied. In the first model the rules are: 

(!} + A + A with probability ti2 

A + A ^ + A with probability t24 

A + id ^ A + A with probability 

A + A ^ A + (l} with probability t^^ 

A + ^ + with probability tig (38) 

$ + A^ A + (l) with probability tga 

— >■ A at the left boundary with probability a 

A — > at the left boundary with probability 7 

A — > at the right boundary with probability (3 

provided that . ^^.t'^'^^.t — = -. — , . . The rules for the second model are: 

A -I- -)• -I- A with probability 1 

-I- A -)■ -I- with probability ti2 

(i) + A ^ A + A with probability t42 ( 39) 

— > A at the left boundary with probability a 

A — > at the right boundary with probability (3 

without any constraints on the microscopic reaction probabilities. It can be shown 
that the steady-states of these models can be written as matrix product states and 
that the dynamics of an appropriately defined single shock front in these models is 
similar to that of a single random walker under parallel dynamics. The results will 
be published elsewhere; however, the classification of the two-state one-dimensional 
driven-diffusive models with open boundaries and nearest-neighbor interactions under 
parallel updating scheme with aforementioned properties, similar to [7] for the models 
under random sequential updating scheme, remains an open problem. 
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